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ABSTRACT 

The mass accretion rate of transonic spherical accretion flow onto compact 
objects such as black holes is known as the Bondi accretion rate, which is deter- 
mined only by the density and the temperature of gas at the outer boundary. A 
rotating accretion flow has angular momentum, which modifies the flow profile 
from the spherical Bondi flow, and hence its mass accretion rate, but most work 
on disc accretion has taken the mass flux to be a given with the relation between 
that parameter and external conditions left uncertain. Within the framework of 
a slim a disk, we have constructed global solutions of the rotating, viscous hot 
accretion flow in the Paczyhski-Wiita potential and determined its mass accre- 
tion rate as a function of density, temperature, and angular momentum of gas 
at the outer boundary. We find that the low angular momentum flow resembles 
the spherical Bondi flow and its mass accretion rate approaches the Bondi ac- 
cretion rate for the same density and temperature at the outer boundary. The 
high angular momentum flow on the other hand is the conventional hot accretion 
disk with advection, but its mass accretion rate can be significantly smaller than 
the Bondi accretion rate with the same boundary conditions. We also find that 
solutions exist only within a limited range of dimensionless mass accretion rate 
m = M /Mb, where M is the mass accretion rate and Mb the Bondi accretion 
rate: When the temperature at the outer boundary is equal to the virial temper- 
ature, solutions exist only for 0.05 rn < 1 when a = 0.01. We also find that 
the dimensionless mass accretion rate is roughly independent of the radius of 
the outer boundary but inversely proportional to the angular momentum at the 
outer boundary and proportional to the viscosity parameter, rh ^ 9.0 aX~^ when 
0.1 < m < 1, where the dimensionless angular momentum measure A = lout/^B 
is the specific angular momentum of gas at the outer boundary lout in units of 
Ib = GM/cs^out, M the mass of the central black hole, and Cg^out the isothermal 
sound speed at the outer boundary. 

Subject headings: accretion, accretion disks — black hole physics — quasars: gen- 
eral — X-rays: general 
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1. Introduction 

The amount of mass gravitationally accreted to the compact objects such as black holes 
is determined by the conditions of gas around the compact objects. The case when gas has 
a spherically symmetric distribution, a polytropic pressure-density relation, and no angular 
momentum was solved by Bondi (1952), who found that the mass accretion rate for the 
transonic accretion is determined only by the density and the temperature of surrounding 
gas, both assumed to approach constant values far from the accreting objects. This rate, 
known as the Bondi accretion rate, is widely used as 'the mass accretion rate' in a variety 
of accretion problems. As a function of the density Poo and the isothermal sound speed Cg^oo 
of gas at infinity, the Bondi rate for pure hydrogen gas is 



y3/2 p3 ' 

where G is the gravitational constant, M the mass of the central object, and 7 the adiabatic 
index of accreting gas. The constant A(7) is 0.25 for 7 = 5/3 and 1.12 for 7 = 1. If we 
define the Bondi radius as vb = GM/c^^^, then Mb — A7~^/^47rr^pooCs,oo- Gravity starts to 
dominate over gas pressure inside the Bondi radius and the Bondi accretion rate is roughly 
the mass flux of gas inf ailing into a sphere of radius with a velocity equal to the sound 
speed at infinity. 

However, many astrophysical accretion flows are expected to have certain amount of 
angular momentum, and the angular momentum will surely affect the flow properties, in- 
cluding the mass accretion rate. Surprisingly, the rate of mass accretion for these rotating, 
viscous accretion flows for given external boundary conditions, has not been studied yet. 
We might expect that if the rotational support is negligible at the Bondi radius, the effect 
of angular momentum to mass accretion rate would be small. But the accretion rate would 
diminish as the dimensionless measure of rotation A = lout/Is approaches unity, where lout 
is the speciflc angular momentum of gas at the outer boundary and Ib = tbCs^oo is the 
representative angular momentum expected at the Bondi radius. 

Although the accretion problem can be reduced to local one in the limit of thin accretion 
disc where the radial velocity is negligible (Shakura & Sunyaev 1973), accretion flow is 
fundamentally global in the sense that the flow structure is only globally determined. As 
pointed out by Yuan (1999), the flow property can be quahtatively affected by the outer 
boundary conditions. The dependence of the mass accretion rate on boundary conditions 
can only be addressed by constructing global solutions. There have been many studies on 
the global solutions of rotating accretion flow (e.g., Narayan, Kato, & Honma 1997; Chen, 
Abramowicz, & Lasota 1997; Nakamura et al. 1997; Lu, Gu, & Yuan 1999) and the effects 
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of outer boundary conditions on accretion flow (Yuan 1999; Yuan et al. 2000), but most 
have focused on the flow structures and emission properties for a given mass accretion rate, 
and the determination of how that rate is flxed by external circumstances remains to be 
addressed. 

In this work, we focus on the mass accretion rate of the accretion flow, especiaUy its 
dependence on the density, temperature, and the angular momentum of gas at the outer 
boundary. Wc do this by constructing the simplest global, transonic solutions within the 
slim disk formalism (Abramowicz et al. 1988). Slim disk formalism uses vertical integra- 
tion to take into account the finite thickness of the accretion disk, and more importantly 
allows for the radial motion of the accretion flow and the critical point unlike the thin disk 
approximation (Shakura & Sunyaev 1973). 



2. Equations and Methods 
2.1. Equations 

In slim disk approximation, the continuity equation is 

M = -AnrHpVr, (2) 

where Vr is the radial velocity {vr < for inflow) and p the density of the flow, both averaged 
over the height of the disk. The disk scale height is if = Cg/flK, where the isothermal sound 
speed = P/p, P is the total gas pressure, and Qk the Kcplcrian angular velocity for 
Paczyhski-Wiita potential (Paczyhski & Wiita 1980). This treatment allows for general 
relativistic effect in an approximate fashion. 

We assume that the gas is composed of hydrogen only and is fully ionized. Ions and 
electrons are assumed to have the same temperature T, and hence in this simplified treatment 
of a pure hydrogen gas, the gas pressure is P = (np + ne)kT = 2npkT = 2{p/mp)kT, where 
Up and He arc the number densities of protons and electrons, respectively, rUp the proton 
mass, and k the Boltzmann constant 

The radial momentum equation is 

Vr^ + {^l-^y+-^=0, (3) 

dr p dr 

where Q(r) is the angular velocity at radius r. Since we adopt Paczyhski-Wiita potential, 
the Keplerian angular velocity is Jl|:(r) = GMr'^[l — {rsch/r)]~'^ where rsch = 2GM/c'^ is 
the Schwarzschild radius and c the speed of light. 
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The angular momentum equation is given by 



(4) 



where the constant a is the viscosity parameter of Shakura and Sunyaev (1973): the viscous 
stress tensor r^r = —OiP (Abramowicz et al. 1988). Depending on exactly how the viscosity 
description is implemented, the parameter can be from 77 = —2 (Abramowicz et al. 1988), 
rj = —1 (Nakamura et al. 1997), 7] = {r /VtK){dVL/ dr) (Narayan, Kato, & Honma 1997) 
to 7] = dhiVLx / d\n.r (Narayan, Kato, & Honma 1997; Yuan 1999). We choose 77 = —1 
because this choice simplifies the angular momentum equation into an algebraic one and has 
an added convenience of automatically satisfying the no-torque condition at the black hole 
horizon (Abramowicz et al. 1988, Yuan et al. 2000). On the other hand, this implementation 
overestimates the shear stress when there is little or no shear, and can cause difficulties when 
the solution is extended to a very large radius. The integration constant Iq is the specific 
angular momentum accreted by the black hole, to be determined as the eigenvalue during 
the construction of solutions for given boundary conditions. 



Finally, the energy equation is 



pVr 



de ^ pd A 
dr dr p 



Q -Q 



(5) 



where e = (7 — 1)~^P/ p is the internal energy of the gas per unit mass, and the heating 
and cooling functions per unit volume, respectively. The only heating process considered in 
the current treatment is the viscous heating, for which we use the description 



= -CaPr(dQ/dr). 



(6) 



We fix = 1 as in Abramowicz et al. (1988), Nakamura et al. (1997), and Yuan et 
al. (2000). Shghtly different values of ( have been used in other works: for example, 
C = -{r/QK){dn/dr) by Narayan, Kato, & Honma (1997) and C = -(din Q/^/dlnr) by 
Yuan (1999). 

Since we mainly deal with high temperature accretion disk with significant radial ve- 
locity, we use the optically thin, relativistic ion-electron bremsstrahlung as the only cooling 
process. 



q — afrl'meC^npne{S2/3){2/n) 



1/2 



kT 



1/2 



1 + 1.78 



kT 



1.34 



(7) 



where et/ is the fine structure constant, Te the classical electron radius, and rUe the electron 
mass (Svensson 1982 and references therein). 
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Equations ([3]) and can be rearranged into the form 

dvr A dcs B 



8 



dr dr Z)' 



where 



D 



B 



A 




pVrC, 



(10) 



(11) 



(9) 



The zero of the denominator, D = 0, yields the well-known critical (or sonic) point condition 
{vr/csY = 27/(7 + 1) (see e.g., Narayan, Kato, & Honma 1997), modified by viscosity. The 
exact value of Mach number Aicr = \vr/cs\rcr the critical point is given by the root of 
the equation 



for our specific choice of viscosity description. For a <^ 1, A^cr^ ~ 27(7 + 1) ""^[l + «^(7 — 



For spherical accretion, the widely used Bondi accretion rate is solely determined by the 
density and the temperature of accreting gas at infinity. To be more precise however, there 
is no unique solution for a given density and temperature at infinity: the same boundary 
condition does allow subsonic (type I), transonic (type II), and unphysical (type III) solutions 
with different mass accretion rates (Bondi 1952). The additional requirement that the flow be 
transonic, i.e., regular at the critical point, becomes an additional constraint, which uniquely 
determines the solution and its accretion rate (cf. Parker 1963). 

In rotating accretion flow, we expect one more physical quantity to be important: the 
angular momentum of the gas. In the slim disk formalism, the integration constant Iq is an 
eigenvalue, and only a specific value of Iq for given M, M, a, pout, and Tout admits a regular 
transonic solution (Muchotrzeb & Paczyhski 1982, Abramowicz et al. 1988). 



(7 + l)Mj - [27 + 2^2(7 - l)]Mcr^ + a'(7 - 1) = 



(12) 



1)(87- 1)(47^)-^]. 



2.2. Boundary Conditions and Method of Calculation 
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Spherical Bondi accretion has a critical point, which prevents the direct integration 
of the equations from the outer subsonic region to the inner supersonic region. This dif- 
ficulty has been handled in two ways: Either one starts with an arbitrarily chosen critical 
radius, integrates both inward and outward with the help of regularity conditions, and ad- 
justs the critical radius until the outer and inner (if any) boundary conditions are satisfied 
(Bisnovatyi-Kogan & Blinnikov 1980; Quataert & Narayan 2000), or one starts with known 
outer boundary conditions such as density and temperature of gas plus initial choice of radial 
velocity (or equally mass accretion rate), integrates inward toward the critical point, checks 
if the velocity suddenly changes sign or diverges, adjusts the velocity (or the mass accretion 
rate) accordingly until the integrated radial velocity behaves normally even very close to the 
critical radius, extrapolates across the critical point when the adjustment is precise enough, 
and then integrates from there inward (Park 1990). This latter iteration procedure has some 
resemblance to the physical behaviors of the accretion fiows. If the gas at outer boundary 
starts to accrete with too small velocity (or mass accretion rate) than that of the transonic 
solution, the gas flow stays subsonic all the way, never passing the critical point (type I). 
On the other hand, if the gas starts to accrete with too large velocity, the velocity of the 
flow diverges before the critical point and becomes unphysical (type III) . Only when the gas 
starts to accrete with the just right velocity, or equally mass accretion rate, can the flow 
pass through the critical point and become steady transonic flow (type II). 

The dynamics of rotating accretion flows are similar except that we must include the 
added complications due to angular momentum and viscosity. The flow must have the 
right mass accretion rate and angular momentum for given density and temperature of gas 
at the outer boundary to pass through the critical point. To find such regular transonic 
solution, the angular momentum eigenvalue Iq is searched by iteratively integrating the 
equations, either from the outer boundary (e.g., Muchotrzeb & Paczyhski 1982, Abramowicz 
et al. 1988, Nakamura et al. 1997) or from the critical point in and out (Chakrabarti 1996; 
Narayan, Kato, & Honma 1997). Similar procedures are also required to find the subsonic 
accretion solutions or supersonic solutions with standing shocks onto compact objects with 
hard surfaces, such as white dwarfs or neutron stars, that satisfy specific outer and inner 
boundary conditions (Popham & Narayan 1991; Narayan & Medvedev 2003). 

Since we are more interested in the mass accretion rate of the transonic flow for given 
boundary conditions at large radius, we integrate from the outer boundary inward with the 
procedure used for spherical accretion as explained above (Park 1990). The basic parameters, 
the viscosity parameter a and the adiabatic index 7 are chosen flrst. In this work, we 
assume 7 — 5/3 and a = 0.01 unless noted otherwise. We then fix the density Pout a-nd 
the temperature Tout of the gas at the outer boundary Tout- We also choose an integration 
constant Iq. We start to integrate inward by arbitrarily choosing an initial guess for the 
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mass accretion rate M, which, along with other conditions, determines the radial velocity 
Vr from equation ([2]) and then the angular velocity Q at the outer boundary from equation 
(jl]) to enable the initial integration. The initial guess for M produces either a subsonic 
solution or an unphysical, diverging solution, and integration stops. We then adjust M until 
the integration can proceed as close to the critical point as possible with a regular velocity 
profile. The bifurcation between subsonic and diverging solution is quite sharp, and the 
iterations can determine M up to an arbitrary precision. Once the integration reaches close 
enough to the critical point, the flow solutions are extrapolated across the critical point to 
the supersonic region by an appropriate rational function, and the integration ensues from 
therein. The interpolation introduces minimal error because the supersonic part of non- 
transonic flow converges to the transonic flow at smaller radii (see e.g., Bondi 1952 and Das 
2007). 

Since fixing Iq and adjusting M is equivalent to fixing M and adjusting Iq and the 
transonic solutions exist only for a rather limited range of Iq, we fix M and adjust Iq in most 
cases. Moreover, equation (jlj) relates Iq to the angular momentum at the outer boundary, 
lout, and, therefore, adjusting Iq is also equivalent to adjusting lout = ^('^ouO^Lt- Since the 
reference angular momentum against which lout needs to be compared is = fsCsirB) = 
GM/cs^out which in fact is the Keplerian angular momentum at if ^ ^schi we define a 
dimensionless angular momentum measure A = lout/ 1 b- To summarize, we construct a global 
transonic solution for given pout, Tout, and A. 

Since our integration starts from a finite radius rather than from infinity, we redefine 
the Bondi rate as Mb = A'nK'^-^/'^{GMf 

PoutCs out- With our choice of cooling function, the 
whole problem can be put into a dimensionless form, and the result, such as the dimensionless 
mass accretion rate m = M/M, is independent of the black hole mass M if the density is 
scaled in units of 

Po ^ - 1.7 X 10-^ (§-) ' g cm-3, (13) 

where ctt is the Thomson scattering cross section. 



3. Result and Discussion 

3.1. Flow properties 

Figure 1 shows typical solutions with different amount of angular momentum: high 
angular momentum (solid lines), intermediate angular momentum (dotted lines), and low 
angular momentum (dashed lines). By high angular momentum, we mean A ~ 1, by low 
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angular momentum, A <S 1, and by the intermediate angular momentum, between the two 
hmits. The density and the temperature of gas at the outer boundary are the same for 
the three solutions: Tout = 2.75 x lO'^K and pout = 1-5 x lO^^po at Vout = lO^r^c/t- The 
dimensionless angular momenta of gas at the boundary and the mass accretion rate are: 
A = 1.7 and m = 0.05 for high angular momentum flow; A = 0.27 and m = 0.386 for 
intermediate angular momentum flow; A = 0.14 and fn = 0.643 for low angular momentum 
flow. Most of the thermal energy is advected with the flow rather than radiated away, and 
temperature stays close to the virial value. 

The flow with high angular momentum becomes supersonic at small radius Vcr = 2.2rsch 
(solid line in Fig. la) and the angular momentum (solie line in Fig. Ic) is close to or a few 
times lower than the Keplerian value (short- long dashed line in Fig. Ic). The radial velocity 
is signiflcant but still smaller than the near free-fall velocity of spherical Bondi flow. This 
is a typical hot, radiatively inefficient advection-dominated accretion flow (ADAF) that has 
been extensively studied (e.g., Narayan & Yi 1994, 1995; Abramowicz et al. 1995; Narayan, 
Kato, & Honma 1997; Chen, Abramowicz, & Lasota 1997). 

In spherical accretion, the flow passes the critical point at a much larger radius ~ r^, 
and the flow becomes supersonic inside ~ r^. But in rotating viscous accretion flow, a large 
part of the flow is subsonic well inside because the rotation of the flow balances out gravity, 
and the value of depends on the angular momentum of the flow. A low angular momentum 
flow (dashed line in Fig. la) has a critical point at a much larger radius = I7rsch 
compared to the high angular momentum one. It has a larger radial infall velocity at the 
outer boundary as well because of smaller centrifugal force (Fig. lb), and the mass accretion 
rate is higher than that of the high angular momentum flow. These characteristics show 
that this low angular momentum flow has more resemblance to the spherical flow than to 
the disk flow as is expected. This type of flow solution is flrst discovered by Yuan (1999) 
and its dynamical and thermal properties have been comprehensibly studied by Lu, Gu, & 
Yuan (1999) and Yuan et al. (2000). 

The flow with intermediate angular momentum (dotted lines) has a critical point at — 
3.2rsch, ^-nd shows intermediate characteristics between the high and low angular momentum 
flow. The flow properties change from disklike to quasi- spherical as the angular momentum 
of the flow decreases. This is similar to the case of inviscid rotating accretion flow, which 
becomes disklike or quasi-spherical, depending on the angular momentum (Abramowicz & 
Zurek 1981). This transition of accretion flow from disklike to quasi-spherical in terms of 
critical radius position is described in detail by Yuan et al. (2000). 
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3.2. Mass accretion rate 

In spherical accretion, the outer boundary is naturally selected to be outside the critical 
point, where the density and the temperature of gas stay roughly constant. In rotating 
viscous flow, the outer boundary is not naturally associated with the critical point which is 
at much smaller radius. So we consider two choices for the outer boundary radius: one at 
the Bondi radius, rB{Tout), for given Tout and the other at a fixed radius regardless of Tout- 
The former choice is equivalent to setting the outer boundary at a virial temperature if we 
define the virial temperature as Tyir{r) = GMmp/{2kr). 

While original Bondi solutions are expressed in terms of the density and temperature of 
gas at infinity, our solutions start at a finite radius. But since the density and temperature 
of gas vary little from infinity to the Bondi radius in the Bondi solutions, we compare our 
solutions for given density and temperature of gas at finite rout against the Bondi solutions 
with the same density and temperature at infinity. 

We first discuss the case where the outer boundary is set to be the Bondi radius, rout — 
GM/c^g^^ oc T^f for given Tout, so that the solutions can be meaningfully compared with 
the Bondi solutions. The density at the outer boundary is either pout = 1-5 x 10~^po or 
1.5 X lO^^po- Since the whole set of equations can be rescaled in density except the cooling 
function q^, the mass accretion rate then is simply proportional to pout when cooling is not 
important. In such case, the flow proflle depends only on m and not on the speciflc value 
of Pout- However, as the mass accretion rate M approaches 0.1 times the Eddington mass 
accretion rate, MEdd = LEdd/(?i cooling becomes important as in spherical accretion (Park 
1990). In such high M / MEdd case, the existence and the properties of the flow depend on 
the specific value of Pout/Po- 

Figure 2 shows the mass accretion rate of constructed solutions as a function of the 
angular momentum at the outer boundary, in rh versus A plane. Symbols represent the 
solutions for different Tout- circles for Tout = 1.1 x 10^° K {rout — 2.5 x lO'^rsch), triangles for 
Tout = 5.5 X 10^ K {rout = 5.0 x lOV^c/.), squares for Tout = 2.8 x 10^ K {rout = 1-0 x lO^rsch), 
pentagons for Tout = 1.1 x 10^ K {rout = 2.5 x lO^rsch), and hexagons for Tout = 5.5 x 10* K 
{'^out — 5.0 X 10 Tsch)- All these solutions are calculated for pout — 1.5 X 10-%. The 
solutions for Tout = 1-1 x 10'^ K {rout = 2.5 x lO^r^c/i), temperature suitable for the ISM in a 
galactic nucleus, have much larger rout, and the bifurcation between subsonic to unphysical 
branch becomes even sharper. Since the critical mass accretion rate above which the hot, 
optically thin accretion disk does not exist also decreases as rout increases (Abramowicz et 
al. 1995) and the Bondi rate increases as Tout decreases, we choose pout — 1-5 x 10~^po so 
that M /MEdd is in a range where solutions can be found. The mass accretion rates of these 
solutions are shown as crosses in the center of Figure 2. 
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We find that regardless of Tout, tlie mass accretion rate m decreases as tlie angular 
momentum A increases. It approaches the Bondi accretion rate as A decreases. This change 
of the mass accretion rate on the angular momentum is expected because given density and 
temperature at Tout, P and H are fixed in the continuity equation (Eq. [2]), and the mass 
accretion rate is determined solely by the radial velocity Vr- From equation Q, the radial 
velocity is determined by the difference between the gravity and the centrifugal accleration 
which is proportional to Z^/r^ (at r ^ tsch)- Larger angular momentum causes smaller infall 
velocity, and the mass accretion rate decreases, and vice versa. 

The mass accretion rate for the lowest angular momentum flow is quite close to the 
Bondi rate (m ~ 1) while that for the highest angular momentum flow is roughly 20 times 
smaller than the corresponding Bondi rate (m ~ 0.05). We also find that m as a function of 
A is rather insensitive to Tout, or equally Tout for m > 0.1: The slope of m(A) has a tendency 
to become slightly flatter (from circles to crosses) as Tout decreases but the difference is not 
large. 

Since we expect the flow profile, especially the radial velocity, to depend on the specific 
value of a, we also calculate the solutions for different values of a. The mass accretion rates 
of a = 0.003 are shown as a series of stars in the left side of Figure 2 and that of a = 0.03 are 
shown in the right side of of Figure 2. Compared to the solutions for a = 0.01 in the center 
denoted by crosses, lower a flows have smaller mass accretion rate and higher a flows larger 
mass accretion rate. Larger a means stronger viscosity for given density and temperature, 
which causes larger radial velocity, and hence larger mass accretion rate, and vice versa for 
smaller a. The relation between m and A for 0.1 < m < 1 can be approximated by 

^(A) = 0.09 (^) X~\ (14) 

The three straight lines (from left to right) in Figure 2 show this fitting function for a = 0.003, 
0.01, and 0.03, respectively. 

The mass accretion rate of the self-similar ADAF that has density pout at rout has the 
same dependency on pout and rout as the Bondi rate Mb (Narayan & Yi 1994). Hence, the 
dimensionless mass accretion rate of self-similar ADAF is simply mADAF = 3.3a. Since 
ADAF generally refers to A ~ 1 flow and m(A = 1) ^ 3m^£)AF, we find that our solutions 
have approximately three times higher mass accretion rates than those of self-similar ADAFs. 
The difference is probably due to the fact that the density and velocity of global ADAFs do 
not exactly follow the self-similar form, especially near the outer boundary (Narayan, Kato, 
& Honma 1997; Chen, Abramowicz, & Lasota 1997). 

For given Tout, the lower limit on the mass accretion rate is determined by the Keplerian 
angular momentum barrier. A smaller mass accretion rate demands smaller radial infall 
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velocity. From equation (jl]), smaller \vr\ requires larger Q, and eventually Q becomes larger 
than the Keplerian value Qk at Vout for too small mass accretion rate. Gas cannot accrete 
in steady-state when the angular momentum is much larger than the Keplerian one, and 
therefore no solution exists when A > 3. Although we could determine m below rh ~ 0.1 
from shooting until A approaches ~ 1.5 — 3 or 3, in many cases we failed to construct the full 
transonic solutions down to the inner boundary. The difficulty is much more severe for small 
a: the bifurcation between subsonic and unphysical occurs at a large radius and computing 
accuracy is not enough to integrate all the way down to the horizon. So we have limited 
confidence in the solutions below m < 0.2 for a = 0.003 and 0.03. The upper limit on the 
mass accretion rate on the other hand is determined by the Bondi accretion rate itself. The 
lowest angular momentum flow has its mass accretion rate already quite close to the Bondi 
rate, which is the mass accretion rate for zero angular momentum flow. Thus any higher 
mass accretion flow will end up on the unphysical branch of solutions in spherical Bondi flow 
(type III), with minimum effect due to the angular momentum. 

We can think of extending the velocity versus radius diagram of Bondi (1952) to a similar 
one with the angular momentum effect added. What angular momentum does is to lower 
the critical mass accretion rate for the transonic flow. As the angular momentum of the flow 
increases, the critical mass accretion rate decreases, and the rotating accretion flow accretes 
transonically at a critical mass accretion rate lower than the original Bondi accretion rate. 
On the other hand, if the angular momentum of the flow decreases, the critical mass accretion 
rate increases up to the Bondi accretion rate which is the critical mass accretion rate for 
zero angular momentum flow, above which no transonic steady-state rotating accretion flow 
exists. Since equation (fT^ gives m > 1 for X < 9a and m cannot be greater than unity, it 
looks like there is no solution for A < 9a. However, we are using slim disk formulation with 
a simple viscosity prescription to describe what is really a two-dimensional quasi-spherical 
flow, especially for lowest angular momentum flow; so the non-existence of solutions below 
the lower limit on A is probably caused by the approximate physical descriptions of the flow, 
especially the overestimated viscosity for the low or no angular momentum flow. Surely, we 
will have a spherical Bondi flow when the angular momentum is zero. Therefore, we expect 
that steady flow probably exists even below A ~ 9a; and 

m ~ 1 for A < 9a. (15) 

We conclude that the steady-state, transonic hot accretion flow exists only when the mass 
accretion rate is within a certain range or the angular momentum at the outer boundary is 
below a certain maximum value. For our speciflc choice of parameters and conditions, the 
range in the mass accretion rate is 

rricr <m< 1, (16) 
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where rricr — Thi^cr) is the mass accretion rate of the flow with maximum possible angular 
momentum at the outer boundary. Our calculations suggest Xcr ^ 3 and rhcr ~ 0.05 for 
a = 0.01. 

Now we discuss flows with the outer boundary temperature different from the virial 
temperature. For example, when the cool, outer geometrically thin disk is attached to the 
hot inner accretion disk, the temperature at the outer boundary will be much lower than 
the virial temperature (Narayan, Kate, & Honma 1997). So we fix Tout = 10"^ rsch but vary 
Tout- The outer boundary density is the same at pout = 1.5 x lO^^po- Figure 3 shows the 
mass accretion rate again in m versus A plane. Symbols are the same as in Figure 2, except 
circles for Tout = 3.6 x 10^ K, triangles for Tout = 2.2 x 10^ K, squares for Tout = 1.1 x 10^ K, 
pentagons for Tout = 5.5 x 10^ K, and hexagons for Tout = 2.0 x 10^ K. 

Again, solutions exist only in a limited range of m, or equally A. But this time both 
upper and lower limits on m vary with Tout- The reason for the difference is that rout is 
fixed and not related to Tout whereas it is related to Tout in Bondi solutions. As in previous 
choice of boundary, the lower limit on rh is from the Keplerian angular momentum barrier. 
However, the upper limit is determined for a different reason in this case. As the mass 
accretion rate for given boundary conditions increases, the radial infall velocity at the outer 
boundary should increase. But the increasing radial infall velocity eventually becomes larger 
than the sound speed at the boundary. The whole flow becomes supersonic from the outer 
boundary to the innermost radius. Since wc arc looking for transonic solutions that uniquely 
determine the mass accretion rate, wc disregard these supersonic branch of solutions, and no 
transonic solutions exist above certain upper limit on m. This upper limit on rh decreases as 
Tout decreases because the sound speed decreases with Tout and the flow becomes supersonic 
at smaller radial velocity, or mass accretion rate, causing rapid decrease of the upper limit 
on rh with decreasing Tout- 

The dependence of m on A is also different. Again, m decreases as A increases, but is 
not inversely proportional to A, nor is independent of Tout- Cooler flow, i.e., lower Tout, has 
lower mass accretion rate m for a given A. When A = 0.1, m for Tout = 2.2xlO^Kis~l 
while that for Tout = 2.0 x 10^ K is as low as ~ 10"^. This means that if gas is injected into a 
given radius with near-Keplerian angular momentum but with temperature much below the 
virial temperature at that radius, the mass accretion rate can be a few orders of magnitude 
below the Bondi accretion rate. 
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4. Summciry 

We have constructed global solutions of the rotating, viscous accretion flow within the 
framework of height-averaged, slim a disk with varying amount of specific angular momen- 
tum A at the outer boundary. We find that the low angular momentum fiow, the fiow with 
A <^ 1, resembles the spherical Bondi flow and its mass accretion rate approaches the Bondi 
accretion rate for the same density and temperature at the outer boundary. The high angular 
momentum flow, the flow with A ~ 1, on the other hand is the conventional hot accretion 
disk with advection, but its mass accretion rate can be significantly smaller than the Bondi 
accretion rate for the same density and temperature at the outer boundary. 

We also find that solutions exist only within a limited range of the dimensionless mass 
accretion rate m, the mass accretion rate in units of the Bondi accretion rate. When the 
temperature at the outer boundary is equal to the virial temperature at that radius, solutions 
exist only when m^,. ^ m < 1 because below m^r the angular momentum of the fiow becomes 
higher than the Keplerian vahie, and the steady accretion is not possible. The exact value 
of TTicr as a function of the viscosity parameter a is not known, but is ~ 0.05 for a = 0.01. 
The upper limit on the mass accretion rate is the Bondi accretion rate since it is the mass 
accretion rate of zero angular momentum fiow. Moreover, we also find that the dimensionless 
mass accretion rate is roughly independent of the radius of the outer boundary but inversely 
proportional to the angular momentum at the outer boundary in units of the Keplerian 
angular momentum at that same radius and proportional to the viscosity parameter, rh ~ 
9.0 aA"^ when 0.1 ^ ^ri < 1. When the temperature at the outer boundary is much less than 
the virial temperature, the mass accretion rate can be a few orders of magnitude smaller 
than the Bondi accretion rate. 

The fact that the mass accretion rate can be smaller than the Bondi accretion rate 
and depends on the angular momentum of gas at the outer boundary could have interesting 
implications for many astrophysical accretion systems, especially underluminous AGNs and 
galactic centers. 

Although we now have a better understanding of the behavior of rotating, viscous 
accretion fiows with varying amount of angular momentum, the precise numerical values 
presented in this work are expected to contain some uncertainties because these hot accretion 
fiows are certain to have real two-dimensional structures, and the one-dimensional, height- 
averaged approximation adopted in this work will only approximately describe the real fiow. 
Hence, more precise analyses of these accretion fiow including other relevant physics such as 
shocks (Chakrabarti & Das 2004), outfiows (Stone, Pringle & Begelman 1999), or radiative 
heating (Park & Ostriker 1999, 2001, 2007; Yuan, Xie & Ostriker 2009; Proga, Ostriker, & 
Kurosawa 2008) will require careful multi-dimensional studies (see e.g. Proga 2007, Kurosawa 
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& Proga 2008, and Proga, Ostriker, & Kurosawa 2008 for current multi-dimensional studies). 
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Fig. 1. — The (a) Mach number, (b) velocity, (c) angular momentum, and (d) temperature 
profiles of typical high (solid lines), intermediate (dotted lines), and low angular momentum 
(dashed lines) solutions. Long-short dashed line in (c) is the Kcplerian angular momentum 
for Paczyhski-Wiita potential. See text for details. 
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Fig. 2. — The mass accretion rate in units of the Bondi rate as a function of the angular 
momentum at the outer boundary in units of Ib for an outer boundary at the Bondi radius 
when the temperature at the outer boundary is equal to the virial temperature. Different 
symbols represent different temperature at the outer boundary: circles for Tout = 1.1 x 10^° K, 
triangles for Tout = 5.5 x 10^ K, squares for Tout = 2.8 x 10^ K, pentagons for Tout = 1.2 x 10^ K, 
hexagons for Tout = 5.5 x 10^ K, and crosses for Tout = 1-1 x 10''' K, but all for a = 0.01. A 
series of star symbols on the left represent Tout — 1-1 x 10^ K and a — 0.003 and those on 
the right Tout — 1-1 x 10^ K and a — 0.03, respectively. 
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Fig. 3. — The mass accretion rate in units of the Bondi rate as a function of the angular 
momentum at the outer boundary in units of Ib for a fixed outer boundary at rout — lO^r^c/i 
and a = 0.01. Symbols represent the same as in Figure 2, except circles for Tout = 3.6 x 10^ K, 
triangles for Tout = 2.2 x 10^ K, squares for Tout = 1-1 X 10^ K, pentagons for Tout — 5.5 x 10^ K, 
and hexagons for Tout = 2.0 x 10^ K. 



